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Abstract 

We consider the effective theories governing the sensitivity to the plasma of certain 
high-energy observables in supersymmetric plasmas, and point out that they preserve 
supersymmetry. Our findings generalize previous observations on asymptotic thermal 
masses in weakly coupled plasmas, to both the real and imaginary parts of self-energies, 
on the light cone and away from it, in weakly and strongly interacting theories. All 
observed supersymmetry violations due to thermal effects turn out to vanish faster 
than E~ 2 in the high energy limit. 



1 Introduction 



Supersymmetry is usually considered as being broken at finite temperature (see e.g. 
PP). A simple reason for this is the different statistics and population functions assumed 
by Bose and Fermi fields, making it hard to see how a Bose-Fermi symmetry could be 
preserved. In Euclidean space, Bose-Fermi symmetry is explicitly broken by boundary 
conditions along the periodic time direction (with period 1/T), which are respectively 
periodic for bosonic fields and anti-periodic for fermionic fields. Nevertheless, one can 
still ask whether the supersymmetry of the underlying equations of motion leaves any 
trace in physical observables. 

One implication of supersymmetry at finite temperature along these lines was de- 
scribed in [2]: due to the existence of a conserved supercurrent, the effective hydrody- 
namics theory which describes the long-wavelength modes of the plasma must contain 
fermionic degrees of freedom. These fermionic modes enter the discussion of non-linear 
effects within this effective theory (e.g. loop corrections), but, of course, are not allowed 
to take on classical expectation values. 

In this letter we propose to look at a different sector of the theory: that of high-energy 
observables. In the strict high-energy limit the plasma decouples and supersymmetry 
is recovered (provided it is present in the vacuum theory). Nontrivial results may be 
obtained by looking at the leading thermal corrections to this limit, such as thermal- 
induced dispersion relations and decay rates. We see no obvious reason why these 
should preserve supersymmetry. Nevertheless, in this letter we wish to report an in- 
triguing fact: for a wide class of high-energy observables, supersymmetry is preserved. 

Our original motivation for this work was to investigate whether the well-known 
observation of supersymmetry for asymptotic thermal masses, in weakly coupled plas- 
mas, did extend to other quantities. We will answer this question in the affirmative, 
and will in fact extend it to all high-energy correlators we could study. More precisely, 
if high-energy supersymmetry violations are characterized by the power of the energy 
E~ n by which they are suppressed (relative to the vacuum correlators), our finding is 
that n > 2 with strict inequality in all considered cases. This is a nonempty statement, 
since the leading thermal effects all have n < 2. 

We will discuss in turn the effective theories we have considered. These include: the 
effective theory for particle masses at weak coupling, in section [2] the effective theories 
for the imaginary part of self-energies (including collinear bremsstrahlung processes 
and 2^2 collisions), at weak coupling, in sectional the self-energies of uncharged 
particles in strongly interacting plasmas with a gravity dual, in section HJ and finally, 
the operator product expansion for deeply virtual correlators, in section O 

We use the phrase "effective theory" to emphasize that the details of the plasma are 
only probed through a restricted set of low-energy correlators, which provide param- 
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eters for medium-independent high-energy effective theories. Supersymmetry should 
be understood as an intrinsic property of these effective theories — we do not believe 
that the thermal nature of the underlying medium plays any role. 

2 Thermal masses at weak coupling 

Thermal dispersion relations (of massless particles) are known to approach the form 
E 2 = p 2 + m 2 ^ at large momenta p ^> gT [3], at the leading order in perturbation 
theory. In applications to supersymmetric theories, it has been repeatedly observed 
that the asymptotic masses moo are the same among particles within a supersymmetry 
multiplet0. Compiling results from the literature [I] [TJ, or by directly evaluating one- 
loop diagrams such as those shown in fig. Q], one arrives at the formulae: 

g 2 C A (Z g + Z)) + g 2 N mattw T M (zf + Z s ) , 
g 2 C M (z g + Z) + Z* + Z s ) + y 2 (z+ + Z s ) , (1) 

where the Zj are certain (non-local) dimension-two condensates that we give shortly. 
Terms in eq. ([TJ are in one-to-one correspondence with particles and interactions of 
renormalizable supersymmetric gauge theories. Ca, Cm and Tm are quadratic Casimirs 
and Dynkin indices for the adjoint and matter representations, respectively, and iV matter 
is the number of chiral superfields, with g the ordinary Yang-Mills coupling. For 
simplicity the Yukawa contribution in eq. ([TJ is normalized to correspond to a term 
~ -^4>ipip + c.c. in the Lagrangian of a single-field Wess-Zumino model. We expect 
supersymmetry to be preserved for more general (e.g. nonrenormalizable) superpoten- 
tials, though we have not checked this explicitly. Nonzero expectation values for the 
D or F auxiliary fields, not considered in eq. (0Q), could break the supersymmetry by 
lifting the bosonic masses^. However, we do not view such effects as being specifically 
thermal sources of supersymmetry violations, since they would do the same thing in 
vacuum. 

The condensates in eq. (JT]), each normalized to give the contribution from two degrees 
of freedom, admit the following gauge-invariant definitions and thermal expectation 

1 Although this observation has been described to me as "well known" in the course of several 
private conversations, I did not succeed in finding a reference to it in the early literature. A recent 
appearance is in [6] . 

2 Is is perfectly consistent to turn on a temperature without generating expectation values for these 
fields. For instance, at the leading order in perturbation theory, F ~ cfycfy vanishes (it involves two 
antiholomorphic fields), and so does D a ~ (f)H a (j) for t a traceless. 
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Figure 1: One- loop fermion self-energy of a fermion if) due to gauge interaction, at 
large energy E. The plasma is probed by soft propagators (with all components ~ T 
in Minkowski spacetime), which give rise to the low-energy correlators ([2]) and are 
denoted here by the crosses. The (hard) remainder of the diagram is expanded in 
powers of T/E. 
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with v M = (1, v) the four- velocity of the hard particle, Ub,f the standard Bose-Einstein 
and Fermi-Dirac distribution functions, and 6a,m the dimensions of the adjoint and 
matter representations, respectively. Useful examples include the thermal masses in 
J\f = 4 SYM, which are all equal to m 2 ^ = g 2 N c T 2 , and the gluon and gluino masses 
in pure glue SQCD, m^ (A) = \g 2 N c T 2 . 

The structures in eq. ([2]) are identical to those entering the hard thermal loop effective 
action [91 ITU] : these are in fact the unique dimension-two gauge-invariant operators 
that can be built out of a light-like four- vector v^ 1 [9]. Given that next-to-leading 
order (O(g)) corrections are associated with hard thermal loop physics [8] (gT scale 
physics), but only ~ g 2 corrections arise from the hard scale ~ E, one concludes that 
at O(g) only the matrix elements <^ get corrected but not their coefficients p]). Thus, 
O(g) corrections preserve supersymmetry. For completeness, we record their values 
here. They are most readily obtained by working within the dimensionally reduced 
Euclidean theory [11], in which the corrections to Z g only arise from the contribution 
Z g — > jj^2A\A\ of longitudinal gluons with zero Matsubara frequency. Next-to-leading 
order thermal masses were also obtained in [12] by means of an indirect thermodynamic 
argument, by relating them to the well-known ~ g 3 T 3 corrections to the QCD entrop)!^. 



3 The authors of p~2] also expressed worries that actual thermal dispersion relations of hard particles 
could grow like m 2 ~ ET at high energies, which would spoil the physical interpretation of the 
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Either way one obtains (Zf LO = Z)° + 0(g 2 T 2 )): 

Z NL0 = II _ lH^L + 0(/T 2 } Z NLQ = ^ _ ^1 + Q{g , T2) (3) 

7rV2 Z7T 

What the factorization formula <^ (and the fact that it preserves supersymmetry) 
becomes when 0(g 2 ) effects are accounted for, if it retains any meaning at all, is 
unknown. The derivation outlined in fig. [1] shows that the supersymmetry of the 
thermal masses at leading order may be interpreted as a statement about the couplings 
of soft particles with various spins to hard propagators: these turn out to be universal 
in super symmetric theories. 



3 Imaginary parts of self-energies at weak coupling 

The imaginary parts of self-energies at weak coupling arise from 2 — > 2 scatter- 
ing against plasma particles, as well as from induced collinear splitting processes 
(bremsstrahlung or pair production). For charged particles in gauge theories, the dom- 
inant contribution to Imll is ~ g 2 TE due to elastic small angle Coulomb scattering, 
though the dominant inelastic contribution arring logarithms) is due 

to induced collinear processes, which we will discuss first. These processes also dom- 
inate the self-energies of neutral particles in gauge theories, provided these particles 
are allowed to split into charged ones. In non-gauge theories, self-energies begin at 
~ g A T 2 E° due to ordinary 2 — > 2 scattering, which we will discuss in subsection 13.21 



3.1 Collinear pair production 

The key aspects of collinear pair production may be briefly summarized as follows. 
This process is only relevant in gauge theories, in which it is initiated by the very 
frequent small-angle (Coulomb) scatterings suffered by either the parent or the daugh- 
ter particles. At high-energies E ^> gT, its long formation time (associated with its 

"effective" thermodynamic masses they have computed. This was motivated by a partial evaluation 
of the gluon self-energy which included only its longitudinal contribution. I believe that such a 
growth is not realized in a complete calculation. Indeed, in the effective theory language we are using, 
the worrisome contribution (A17) in [12j can be matched to a matrix element of the dimension-1 
gauge-invariant gluon condensate iv-A8{— iv-D)v- A, which arises from the eikonal coupling of charged 
particles to gauge fields. Although this operator definitely contributes to the imaginary part of self- 
energies — it gives rise to the universal T ~ g 2 T log(l/g) damping rate of charged particles [TJ] — as I 
will discuss more fully in a forthcoming paper [13j , the longitudinal and transverse gluon contributions 
exactly cancel out from its real part. This cancellation can be foreseen through a causality argument: 
the real part (of a time-ordered correlator, which is the natural operator ordering at high energies) is 
equal to that of a retarded correlator, which cannot probe the plasma at light-like separations. Thus, 
the thermodynamic results of [12] ought to agree with eq. ([3]). 
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collinearity) requires that multiple thermal scatterings occurring during it be summed 
coherently. This causes a parametrically significant destructive interference, the so- 
called LPM effect [15j . which is responsible for the non-analytic behavior II oc E 1 ! 2 . 
For relativistic plasmas, a complete leading-order treatment was given (for photons) in 
[T6] (see also earlier discussions [T8l fT7] . in which different approximations are made). 
Somewhat schematically, the result may be written in the form: 

- 2ImII a (£) = V / dzP a ^ hc {z)F a ^ bc {E, z) (4) 
be Jo 

where P a ^bc are ordinary DGLAP kernels [20J, governing collinear physics, be indexes 
final states and z = Eb/E a is the longitudinal momentum fraction. In eq. (jl]) we have 
omitted final state Bose-enhancement or Pauli-blocking factors which are not needed 
unless z or (1 — z) are very small, ~ T/E. The functions F(E, z) depend in a compli- 
cated way on E and z and are to be obtained by solving an effective inhomogeneous 
Schrodinger equation governing the evolution of the pair in the transverse plane |16j . 
This equation depends on the details of the plasma through a collision kernel C(q±), 
which is a function of the transverse momentum transfer. Its only property that we 
need is that it involves only eikonal physics: it does not depend on the spins of the 
particles. For our purposes F(E, z) in eq. PJ is thus just some universal function, 
the same for all final states among a given supersymmetry multiplet. In the leading 
logarithmic approximation [IE] . F(E,z) ~ g A N^T^ E^ z^^ (1 — z)~* (l g(^^— ^y)) 1 / 2 . 

The only ingredients in eq. (jlj) which could possibly break supersymmetry are the 
DGLAP splitting kernels P a _ + b c (z). Such kernels are listed in table [U for various 
supermultiplets of initial and final states. As shown in the table, when complete 
supermultiplets of final states are summed over (thereby enforcing the symmetry under 
z — > (1 — z)), supersymmetry with respect to the initial particle is restored. Not 
shown in the table (it is related to the first three entries by a crossing symmetry 
|21j). but which also preserves supersymmetry, is the process of bremsstrahlung of 
a gauge particle off a chiral multiplet. Thus, all in- medium splitting rates preserve 
sup er sy mmet ry. 

Observations of supersymmetry in DGLAP kernels were made long ago by Dokshitzer 
[2T] . and subsequently given an explanation (at the one- loop level, we believe) by 
Lipatov and collaborators [22]. Here we are merely reporting their implications at 
finite temperature. 

We expect coupling constant corrections to eq. (jlj) to first arise at O(g). In thermal 
perturbation theory such ~ g factors arise from ordinary loop factors g 2 multiplied by 
large bosonic occupation numbers % ~ T/p° ~ T/gT, and are associated with hard 
hard thermal loop physics [8] at the gT scale. Such soft physics can only interfere with 
processes having a sufficiently long duration, such as the soft scatterings contributing to 
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Process 


DGLAP kernel P(z) 


Sum 


7 — > ip^ijj 
7 — > (jy(p 
7 — > (jyij) 


e 2 [2 2 + (1 - zf] 
e 2 [2z(l - z)] 


e 2 


e 2 [2z] 


e 2 


9 ^ 99 
g -> At A 
A^<?A 


g 2 C A [2(1 - z)/z + 2z/(l -z) + 2z(l - z)\ 

g 2 CA [z 2 + (1 _ z) 2] 


<? 2 C A [2/^ + 2/(1-^-3] 


<fC A [4z/(l-z)+2(l-z)] 


(7 2 C A [2/^ + 2/(l-^)-3] 


— >• 

•0 — > 0t-0T 


?/ 2 [l] 




2/" [2*] 


2/ 2 



Table 1: DGLAP splitting kernels. The first three entries govern the splitting of 
photons and photinos to a chiral multiplet, the next three pertain to the Yang- Mills 
splitting of a nonabelian gauge multiplet, and the last two govern the Yukawa splitting 
of a chiral multiplet. Supersymmetry is restored when complete supermultiplets of 
final states are summed over. 

C(q±) with qj_ ~ gT, so we believe that this is the only ingredient suffering from 0{g) 
corrections. These soft collisions have a purely diffusive effect, so we could equivalently 
say that all 0(g) corrections pertain to the so-called transverse momentum diffusion 
coefficient u q" . These corrections should be calculable using techniques similar to those 
used for heavy quark diffusion in [23], but this has not yet been done. Since only eikonal 
physics is involved, these corrections trivially preserve supersymmetry. 

The 0(g 2 ) corrections to eq. (@J are expected to possess a much more interesting 
and richer structure. For instance, they will most certainly require dealing with 
the scale dependence of the partonic constituents of the plasma, which should ul- 
timately lead to "saturation" effects [21] at very high energies upon summation of 
large logarithms a s log(E/T) and a s log(q 2 _/T 2 ) with q\ ~ E 1 ! 2 ^! 2 . The scale evo- 
lution of the constituents of the probe, which has to be treated in the presence of 
the LPM effect, should also enter at this order. Other interesting (though manifestly 
supersymmetry-preserving) effects may include sensitivity to nonperturbative g 2 T-scale 
magnetic physics, which we believe contributes to q at 0(g 2 ). We leave to future work 
a detailed analysis of these effects and of the question of whether they preserve super- 
symmetry. 

As for the subleading corrections in T/E, we expect supersymmetry-breaking effects 
in II not to be larger than ~ T b l 2 E~ l l 2 . These could arise from various ~ T/E or 
~ q\j 1 E 2 ~ (T/E) 3 / 2 corrections to ingredients entering F(E,z), such as the eikonal 
vertices. 
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Process 


|M| 2 /V 


Processes 


|M| 2 /V 


ijjijj — > ifjijj 

00 -> 00 
0-0 — > 0-0 


1 
1 

—u/s 


ipip -> 0V> -»• x 

^ -> X, 0^ -> X 
^0 -> X, 00 -> X 
$<j>->X,<j><j>->X 


1 

[2+f+a 
i 

[2 + f + ^] 



Table 2: Left panel: scattering amplitudes \Ai\ 2 in Wess-Zumino model, with am- 
plitudes related by crossing symmetry not shown. Right panel: amplitudes summed 
over final states, for which supersymmetry is manifest as a function of particle 1 with 
particle 2 held fixed. 



3.2 2^2 scattering at weak coupling 

Ordinary 2-^2 collisions dominate self-energies in non-gauge models, which we will 
now discuss; their total rate is found to preserve supersymmetry. We first recall the 
general formula for the total collision rate (— Imll = — TE): 

x \M 1S2 ^ S3S4 \ 2 n b (E 2 )(l±n c (E 3 )){l±n d (E 4 )). (6) 

S2S3S4 

Here the particle labels are as defined as in fig. (2J the Sj label the corresponding particle 
species, and rii are the corresponding distribution functions. 

Let us first assume, for a moment, that the distribution functions can be omitted 
in the final state ( "Bose-enhancement" and "Pauli-blocking" ) factors (1 ±rij), which 
is justified for generic final state energies E3 ~ E4 ~ \fE\E~i ~ yET. The integrand 
then depends only on the sum ^ S3S4 \M\1 S2 _, S3Si . Such matrix elements summed over 
final states turn out to obey supersymmetry identities, with respect to the particle 1, 
for fixed identities of particle 2. This is exemplified in table [2] for single- field Wess- 
Zumino model with cubic superpotential, and the generalization to other models will 
be given shortly. Therefore, the contribution to eq. ([6]) from the region E 3l E^ 3> T 
preserves supersymmetry. 

It is easy to convince oneself that for bounded amplitudes |.M| 2 , the regions E 3 ~ T 
and E4 ~ T suffer from ~ T/E phase-space suppressions, and so the neglect of the 
final state distributions in eq. (jSJ) is justified. However, s/t ~ ET/T 2 singularities in 
squared matrix elements when t < T 2 can overcome this suppression, and a separate 
discussion is required for these singular termfl (The physical equivalent region u — > 0, 



4 The total integral of such ~ l/t singularities is logarithmically divergent at t — ► 0. This is cured 
by resumming hard thermal loop self-energies [8] to the soft exchanged fermion propagator. 
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Figure 2: 2 — > 2 scattering processes in Wess-Zumino model; solid lines are fermions 
and dashed lines are scalars. 

or E 3 ~ T, can be treated similarly). To establish the supersymmetry of this region, 
for which the distribution function n{E^) must be kept, we need another ingredient: 
universality of the 1/t singularities. Indeed, the coefficient of 1/t at t — *■ 0, which is 
due to soft fermion exchange, is left unchanged when the hard particle 1 is replaced by 
its superpartner (e.g. if particles 1 and 3 are exchanged in fig. [2]). This shows that the 
complete ~ T 2 E° self-energies in the Wess-Zumino model preserve supersymmetry, up 
to ~ T^E^ 1 corrections. 

It is interesting to analyze the t — > singularities by means of the effective theory 
language used for thermal masses in section [21 Indeed, the region E^ ~ T 1 t ~ T 
in fig. [2] is characterized by soft fields coupled to a hard line and is thus governed by 
the gradient expansion depicted in fig. [TJ This means that the ~ T 2 E° contribution 
to eq. ([6]) from soft fermion exchange may be understood as a contribution to the 
imaginary part of the fermion condensate in eq. (T5]), at one-loop in thermal perturbation 
theory^. This way we see that the universality of 1/t singularities as a function of 
particle 1 is closely linked to the supersymmetry of thermal masses. 

We now show, as promised, that supersymmetry of scattering amplitudes summed 
over final states holds in any supersymmetric theory, as a property of the S-matrix. 
Introducing the notation Pi 1 ...i n for projection operators which perform the sum over 
complete supermultiplets of scattering states with n particles, this follows from con- 
sidering the following trace (over scattering states): 

Tr[tfP u S(\2)(2\® [Q,|1)<1|])], (7) 

with 5* the S-matrix and 1 denotes the superpartner of particle 1. For any supersym- 
metry generator Q which does not annihilate particle 1, the commutator [Q, |1)(1|] oc 
(|1}(1| — |1)(1|), so eq. ([7]) computes the difference: 

^(|M 12 _ 3S4 | 2 -|M i2 _ 3S4 | 2 ). (8) 

5 For scalar exchange, this argument relates the locality of the scalar condensate (j>*<j> in eq. ([2]) 
(which implies that it is purely real) to the absence of 1/t singularities due to scalar exchange. 
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For a massless particle 2 one can always choose Q so as to annihilate particle 2; such 
a Q commutes with |2)(2|, with the S-matrix as well as with the projectors Pi x ...% n 
(by construction), showing that eq. (jTJ) (and thus eq. (jH])) vanishes, being the trace of a 
commutator. This establishes the supersymmetry of the contributions from E 3 , E^^> T 
to eq. (jSJ) in any theory. 

Combining the results of the preceding sections, we have reached a simple conclusion: 
the full thermal self-energies of gauge-neutral particles preserve supersymmetry, at 
leading order in the coupling, up to corrections suppressed by a at least T 5//2 i? _1 / 2 . 
Although we believe the analysis can be generalized to charged particles (for which the 
analysis is made more complicated by the stronger singularities Ai ~ 1/t associated 
with gluon exchanged, and various sources of infrared divergences which make these 
self-energies less cleanly defined), here we will refrain from doing so: we are content 
with a robust result for gauge-invariant self-energies. 

4 Strong coupling 

Maldacena's conjectured gauge/gravity correspondence [27] renders possible, among 
other things, the calculation of correlators of currents in certain strongly coupled large 
iV c gauge theories. In theories which have a continuous R-symmetry, such as the 
SU(4) of M = 4 super Yang-Mills, "photons" and "photinos" can be introduced by 
weakly gauging a U(l) subgroup of the R-symmetry. Their self-energies are then given 
by suitable two-point functions of currents and their superpartners, which we now 
evaluate by means of the correspondence. 

In the case of the on-shell photon self-energy in TV = 4 SYM, it was argued by 
means of a WKB approximation [28] (in appendix) that at high energy the calculation 
localizes itself near the boundary of the AdS space. Here we generalize this phenomenon 
to other backgrounds, and establish two properties of the resulting effective theory 
of high-energy photon/photino propagation. First, it only probes the underlying low- 
energy medium through the expectation value of its energy-momentum tensor (actually, 
only through one component oc p^p u T^ v ). Second, it preserves supersymmetry: the 
absorption rates and dispersion relations of a photon and of a photino are identical. 

6 In an effective theory language this implies the existence of certain pure imaginary operators 
which are not probed by the thermal masses, such as the dimension-1 operator described in footnote- 
foot:diml. At dimcnsion-2 I find operators like v-AS(—ivcdotD)D (representing e.g. an interference 
term between t-channel gluon exchange and Z?-term scalar self-interaction in (jxjj — > (jxj) scattering), as 
well as its superpartner involving A. At least for vanishing chemical potential I could check that such 
operators take on no expectation value. 
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We will be considering five-dimensional metrics of the general form 

dg 2 _ R 2 9(z)dz 2 + h IMU (z)dx tJ 'dx u 
z 2 

for which, near the boundary z = 0, the metric approaches that of AdSs with radius 
R (for which g(z) = 1 and h^z) = r/^). The metric eq. (EJ) should be sufficiently 
general to cover any system invariant under space-time translation that admits a gravity 
dual. For the AdSs black hole, relevant for M = 4 SYM at finite temperature T, 
— /loo = 1 _ (kTz) 4 , hij = 5ij, h i0 = and g(z) = (— /loo) -1 - At certain steps below 
rotational invariance will be assumed; these steps will be highlighted. 



4.1 Bulk equations 

The bulk dual of the spin-1 current which couples to the photon is a five-dimensional 
gauge field, whose field strength tensor obeys Maxwell's equations: 



, 0, I h ~, + , (10) 

i/g(z)det(-h(z)) \ z \j g(z) ) 

^aFfjiu 9^F ai , d u F a ^ . (11) 

Here /x, u, a, p are space-time indices but a may cover all five coordinates. We will 
restrict our attention to space-time momentum eigenstates = ip^. A closed equation 
for the transverse electric field F ±, for v =_L a component perpendicular to p M , may 
be obtained by acting on the first equation with a partial time derivative do, and 
using the second equation. Specifically, one uses relations such as d F z ± = d z F ±, 
which follow from dropping perpendicular derivatives d± in the latter. To fully exploit 
such simplifications, rotational invariance must be assumed, so that upstairs derivative 
h ±cr d a also vanish. This yields the closed equation: 



Zk±± -A f ^J'^P^d^ | = W-wvPol, (12) 



^g(z)det(-h(z)) \ z V g(z) 

in which no summation over _L-indices is implied. 

The bulk dual of the spin-| operator coupling to the photino is a five-dimensional 
Dirac fermion with mass m — | [26] (in units with R — 1). It possesses as many 
components as two four- dimensional Weyl spinors, but it is dual to only one such 
spinor: the sign of m breaks the symmetry between the two Weyl components. The 
bulk Dirac equation reads: 



[p + m]ip = 



V>, (13) 
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with a, a = ... 4 and e" the orthogonal basis. Under the assumption of rotational 
invariance, the term involving the spin connection uj must be proportional to the single 
matrix j z , and it can be removed by a z-dependent field rescaling. We choose the 
rescaling tp = z 2 (det(— h))~ 1 / 4 e~ m f dz \T^l z ^)^ which leads to the following equations 
for the Weyl components of ipL,R of i/j: 



2m 
-o z 

z 



R, 



(14) 
(15) 



Here j>L,R are Weyl operators associated with the four-dimensional metric h^ u (z). With 
m = +| the component relevant near the z = boundary is ipL and we are calculating 
the self-energy of a left-handed photino. Eq. (fT5l) implies a closed equation for ip L : 



1 



9{z) 



--d z 



2m 

z 



1 



RV9{Z) 



L ■ 



(16) 



4.2 WKB solution and super symmetry 

We are now in position to discuss the WKB approximation. By means of a change of 
variable y = y(z), Maxwell's equation eq. (1121) may be cast in a Schrodinger form with 
potential proportional to the squared energy p%, provided 



d ^ = 2zh x J A ?W (17) 
dz V det(-h(z)) v ; 

The factor of 2 is convenient: near the boundary y ~ z 2 . For black holes (like the 
AdSs black hole metric given above) the function g(z) has a pole at a finite value of z 
(the location of the horizon), while the function det(— h) vanishes there. In this limit 
y is mapped logarithmically to infinity. The rescaled potential remains finite there, 
though, due to the blowing up of h 00 ; it then depends only on the energy E = p°. 

The qualitative features of the Schrodinger potential in the equation [d 2 — V(y)]F ± = 
are sketched in fig. [31 The shape of the potential depends on the geometry but not 
on the energy, which only determines its overall normalization. At large y the potential 
becomes constant, while for y —>■ the leading term becomes, for on-shell and off-shell 
momenta respectively, 

V (V)->1 y dhr(z) n 2 T^p,p u „ ( 18 ) 

4 P ^^T = - y 2N 2 ' P = °- 
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~l/m 



Figure 3: Schematic features of the Schrodinger potential V(y)/po, when p 2 = 0. It 
approaches the universal linear behavior eq. (1181) near the boundary and tends to a 
constant at the horizon y — > oo, with a transition regime that may depend on the 
details of the theory and on possible intrinsic mass scales m. 

Here we have used that the leading corrections to the metric near the boundary are 
proportional to z A and are related to the expectation value of the stress-energy tensor 
T^; its trace part, if nonzero, does not contribute when p 2 = 0. The normalization in 
eq. (fl8l) is appropriate to the Af = 4 SYM theory. 

At the horizon y — > oo the solutions are oscillatory and in-falling boundary conditions 
Fox e%w must be imposed for calculating retarded correlators [25], with u = p°7rT/2 
the natural frequency near the horizon. To obtain correlators of currents, as described 
shortly, this solution must be continued to the AdS 5 boundary z = 0. For energies high 
relative to all intrinsic scales in the metric, a WKB approximation can be used. This 
is applicable for y down to y ~ 1/p 2 (resp. y ~ (T 4 _E 2 ) -1 / 3 ) for p 2 (resp. p 2 = 0), 
at which it breaks down due to the redshift factors^. Note that these scales are the 
intrinsic scales of the Schrodinger equations with approximate potentials ffl8|) . The 
problem is thus reduced to solving those approximate equations, with large y behavior 
matching the WKB form oc y-^e^" dyVV . 

For the Dirac equation eq. (I16p one needs only modify the change of variable eq. (117)) . 
to: 

? = 2zy/fc)J, R{Z) v (19) 

dz yyy } M^ = o) y } 

For the on-shell component ^£ of a left-handed photino in a rotationally-invariant 
background, the operator fa is nonsingular with eigenvalue E(^/\h 00 \ + \//i 33 ). Here 
/i 33 gives the metric component along the longitudinal direction. Just as for eq. ([Tj 



7 The transition between the two regimes p 2 ^ and p 2 = takes place smoothly around \p 2 1 ~ 
E 2 / 3 T 4 / 3 , for which value of p 2 the two estimates for y cross each other. This natural scale p 2 (E) is 
exactly the "saturation scale" p s ~ T/x s discussed in |29j at strong coupling, viewed as a function of 
E with ie s =p 2 /2ET. 
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near the boundary y ~ z 2 and the horizon is mapped logarithmically to y = oo, so the 
same WKB approximation applies. More significantly, one readily sees from eq. (1161) 
that the approximate potentials near the boundary are identical to the photon case, 
eqs. pi). 

Correlation functions are obtained by prescribing the limiting values of the fields 
F ± and ipi near the boundary and evaluating boundary terms oc d y F ± (see e.g. [28J), 
or proportional to iptp ~ ^r/z ~ -j^dyipL |26J, or more precisely, 

n 7 = -^,i m ^iM, n^J^ta^M ,20) 



■in' 



F 0± (y) ' "" ' 8tt 2 v-k) ip-£(y) 



Here II 7 = uLu is the photino self-energy sandwiched between on-shell polarization 
spinors u, whose real part yields the thermal mass squared. The normalizations eq. (1201) 
are fixed by the well-known supersymmetry-preserving vacuum result, Il 7 = IT 7 = 
—N 2 p 2 /327i 2 log(p 2 /fj, 2 ), p 2 = p 2 — p%. On the light-cone, Schrodinger's equation with 
the approximate potential ( fl8l) is solved in terms of Bessel (Hankel) functions, 

F ±(y) = i>l=y 1 > (j$ (j^ f J + ^ Q^ 1 )) > ( 21 ) 

with u 2 = ir 2 T IMU p > *p u /2N 2 , yielding with eq. $3$ the result: 

n 7 (p) = IT 7 (p) = ig ^ r 1 ^ (3l - i3§) J, p°=p large. (22) 

The imaginary part of this result reproduces that given in [2S], with u = p°T 2 7c 2 /2 in 
thermal M = 4 SYM. Corrections in T/E to this result may be found by expanding 
the potential eq. (|18p to higher orders near the boundary; for the AdSs black hole this 
expansion proceeds in powers of y 2 ~ co>~ 4 / 3 , so the first subleading corrections to II 
are ~ lu~ 2 / 3 . 

It is remarkable that photon self-energies at strong coupling, and high energies, 
depend only on one property of the plasma: its stress-energy tensor. An heuristic 
picture of strongly coupled plasmas, based on the idea of parton saturation, has been 
proposed recently [29] in which such a universality comes out naturally. 



5 Operator Product Expansion 

Supersymmetry violations in deeply virtual correlators may be analyzed by means of 
the operator product expansion (OPE) [30J. The OPE is a means of separating short- 
distance and long-distance physics, so that the thermal corrections to deeply virtual 
(short- distance) correlators with E ^> T are expressed in terms of the expectation 
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value of local operators. Thermal corrections are thus suppressed by powers ~ E A , 
with the powers A set by the scaling dimensions of local operators^. 

The difference between a correlator of operators and of their superpartners may be 
written as a supersymmetry variation (which corresponds to the fact that it vanishes in 
supersymmetry-preserving vacua). For instance, for correlators of transverse currents 
e- J and of their superpartners X a [32J, one schematically has: 

evJip) e T J - l -\\p) i x i>i 2 X oc efQ a [a| (p)e 2 -] (23) 

with = an d a, a spinor indices. As an operator equation the OPE must 

commute with the supersymmetries, so from the OPE of the right-hand side of eq. (1231) 
one concludes that the operators on its left-hand side are supersymmetry variations. 
Thus, in order to see supersymmetry violations at order E~ 2 or stronger, one must find 
local gauge-invariant fermionic operators of dimension | or less. 

In a wide class of theories, an accidental symmetry comes into play: such operator do 
not exist. These theories certainly include all weakly coupled gauge theories containing 
no U(l) vector multiplets and no gauge-neutral chiral superfields. In these theories, 
gauge-invariant dimension-2 bosonic operators (such as Tr00 or Tr0*</>) do not corre- 
spond to any supersymmetry variations, and must thus enter in the same way in the 
OPEs of fields and of their superpartners; this is easily verified at leading order for the 
OPE of currents in weakly coupled gauge theories. The lowest-dimensional fermionic 
operators are dimension- 1 supercurrents, from which we conclude that thermal super- 
symmetry breaking can only be seen through dimension-3 operators, ~ E~ 3 . 

When neutral chiral superfields or U(l) vector multiplets are present, nonzero expec- 
tation values for D ~ cjfcj) or F ~ (jfcff auxiliary fields (which enter the supersymmetry 
transformations of the gaugino and fermionic matter fields, respectively) could produce 
supersymmetry violations at dimension 2. A similar possibility was observed for ther- 
mal masses, in section eq. (T5]) but, as we discussed there, we do not view it as being 
specifically related to thermal effects. Thus, we conclude that in weakly coupled theo- 
ries, the absence of thermal supersymmetry breaking below dimension 3 is generic. 

It is not possible to analyze general theories at finite values of the coupling constants, 
because finite anomalous dimensions can alter the power-counting. Nevertheless, for 
certain strongly coupled theories accessible to the AdS/CFT correspondence, it is easy 
to be more quantitative. For instance, it is known [33J that in Af = 4 SYM at large 't 
Hooft coupling A > 1, only protected (chiral) operators have finite dimensions A <C 
A 1//4 and that the lowest-dimensional fermionic operator has dimension 2+ 1 = | (being 



8 It is worth noting that the OPE, being most rigorously formulated in Euclidean signature, 
does not enjoy at the moment the same rigorous status in Minkowski space-time with respect to 



nonperturbative physics (see e.g. [3T]). 
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the supersymmetry variation of a dimension- 2 primary). Similarly, the M = 1 theory 
dual to IIB string theory on AdSs x T 11 [M] is known to contain no fermionic operator 
of dimension less than 2 [35]. Thus, in these theories, supersymmetry violations can 
only be seen at the level of ~ E~ 3 or ~ corrections, respectively. A discussion of 
more general strongly coupled theories will not be attempted here. 

6 Discussion 

In this letter we have shown that supersymmetry is a generic property of the effective 
theories which describe high-energy correlators in supersymmetric plasmas. These 
correlators included self-energies at high energies on the light-cone as well as far away 
from it (large virtuality). 

For all correlators (except for the more tractable deeply virtual correlators, treated 
in section [5]) our analysis has been limited to the leading nontrivial order in the relevant 
coupling constant expansion (at both weak and strong coupling). Unfortunately, this 
seems like a major limitation of this work: it makes it hard to decide whether our 
findings highlight general properties of supersymmetric theories, or whether they are 
mere artefacts of these extreme limits. In this respect it would be interesting to know 
the structure of higher order corrections, at both weak and strong coupling. 

We have found that thermal supersymmetry violations in all correlators are sup- 
pressed by a power of the energy E~ n (relative to the vacuum correlators), with n 
strictly greater than 2. (Violations with n = 2 were observed in sections [2] and [5] due 
to nonvanishing D-term of F-term expectation values, but we do not regard these as 
specifically thermal effects. These would have identical supersymmetry-breaking effects 
in vacuum correlators.) 

We find pleasing that such a simple uniform bound holds. This makes one wonder 
whether our findings could be elevated to some sort of theorem, valid independently of 
perturbation theory, though we have no concrete proposal to make along these lines. 
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